Transport through a quantum spin Hall quantum dot 
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Quantum spin Hall insulators, recently realized in HgTe/ (Hg,Cd)Te quantum wells, support topo- 
logically protected, linearly dispersing edge states with spin-momentum locking. A local magnetic 
exchange field can open a gap for the edge states. A quantum-dot structure consisting of two such 
magnetic tunneling barriers is proposed and the charge transport through this device is analyzed. 
The effects of a finite bias voltage beyond linear response, of a gate voltage, and of the charging 
energy in the quantum dot are studied within a combination of Green-function and master-equation 
approaches. Among other results, a partial recurrence of non-interacting behavior is found for 
strong interactions, and the possibility of controlling the edge magnetization by a locally applied 
gate voltage is proposed. 
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I. INTRODUCTION 

Topological insulators and superconductors*— have re- 
cently become the topic of extensive experimental and 
theoretical research. These materials can be described 
in terms of weakly interacting quasiparticles and their 
single-particle Hamiltonians show non-trivial topological 
properties in momentum space. They have an energy gap 
in the bulk but a topologically protected gapless spec- 
trum of boundary states. Schnyder et al£& and Kitaev 5 
have put forward an exhaustive classification of these 
systems in terms of their Altland-Zirnbauer symmetry 
classes^ and of the number of spatial dimensions. 

A topologically nontrivial state is possible for the sym- 
plectic class All in two dimensions, corresponding to sys- 
tems with spin-orbit coupling in the absense of an ap- 
plied magnetic field and of superconductivity. This so- 
called quantum spin Hall (QSH) stated has protected 
edge states with gapless Weyl-type dispersion. Bernevig 
et al& have predicted and Konig et aL 10 ^ 11 have observed 
the QSH effect in HgTe/(Hg,Cd)Te quantum wells. 

The edge states of the QSH system show spin-mo- 
mentum locking in the sense that right-moving (left- 
moving) electrons are strictly spin- up (spin-down) ji&II 
Consequently, density-density interactions cannot lead to 
backscattering since they cannot flip the spin. On the 
other hand, spin-dependent scattering, for example by 
magnetic impurities, can lead to backscattering. Uncon- 
ventional transport properties are thus expected and pos- 
sible applications in spintronics can be envisaged. It is 
therefore of interest to study electronic transport in pro- 
totypical device geometries involving QSH edge states, 
which is the subject of this paper. 

The underlying idea is to realize a quantum dot as 
a finite-length segment of a QSH edge. This cannot be 
achieved by electrostatic gating since an electric potential 
just shifts the edge bands up or down without opening 
a gap and thus does not lead to the formation of tun- 
neling barriers. However, such barriers could be realized 
by ferromagnetic insulators grown in contact with the 
edge. They would impose an exchange magnetic field 



FIG. 1: (Color online) Cut-away view of the QSH quantum 
dot. 



orthogonal to the spin-orbit field and open a gap. We 
also include a gate electrode that can be used to tunc 
the electrostatic potential on the quantum dot. Figure Q] 
shows a sketch of the device. Since the QSH edge states 
typically have a small Fermi wavenumber kp, we consider 
the case that the barrier width is negligible compared to 
Xp = 2n/kp. We restrict ourselves to the cases of par- 
allel and antiparallel exchange fields in the two barriers. 
Finally, we assume the thermal energy kpT to be small 
compared to the lifetime broadening of the dot levels so 
that we can set the temperature to zero. Including a fi- 
nite temperature is straightforward. An alternative real- 
ization of a QSH quantum dot as the edge of a small QSH 
puddle has been analyzed by Tkachov and Hankicwicz, 12 
for negligible charging energy. They have also studied 
contributions to the resistance due to the coupling of a 
QSH edge to normal leads and the effect of an orbital 
magnetic field parallel to the spin-orbit field. 12 

After introducing the model in Sec. [H] we study the 
effect of the gate voltage and the bias voltage on the 
transport in Sec. IIII1 A Landauer approach 13 is used to 
obtain the current for arbitrary strength of the magnetic 
barriers. We will see that the result can also be analyzed 
in terms of a generalized Meir-Wingreen (MW) formula 
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for the curren t 14 i 15 in terms of the non-interacting Green 
function. It is generalized in the sense of pertaining to a 
real-space continuum Hamiltonian instead of a tunneling 
Hamiltonian, which leads to a nontrivial form of the dot- 
lead coupling functions T L , T R . 

In the next step, the effect of electron-electron inter- 
action on the dot is studied in Sec. IIV1 We again employ 
a simple model by including a charging energy in terms 
of the excess charge. This is valid if the range of the 
electron-electron interaction is large compared to the dot 
size W, i.e., for small dots, but should give qualitatively 
correct results beyond this regime. The opposite case of 
short-range interactions is certainly of interest. For an 
unmodulated QSH edge this case has been studied for 
example in Refs. [HI, [T3, EH- A quantum dot in a spinless 
Luttinger liquid, not a QSH edge, has been studied by 
several authors.— ~— Tunneling through a quantum dot 
between two QSH edges described as Luttinger liquids 
has been addressed by Law et al£2- 

In Sec. IIV1 we employ an equation-of-motion approach 
to obtain the multiple-pole structure of the spectral func- 
tion appearing in the generalized MW formula. We dis- 
cuss possible approximations for the weights of the poles 
and show results calculated from a master equation in the 
sequential-tunneling approximation with additional level 
broadening, valid for weak coupling through the magnetic 
barriers. We summarize the main results in Sec. [V] 



II. MODEL 

The non-interacting part of the Hamiltonian of the 
QSH edge with magnetic barriers is written as a 2 x 2 
matrix in spin space, 

H = —ihvF<y z d x — rjhvFO- x S(x) =F r]hvFcr x S(x — W) 
+ V(x), (1) 

where vf is the Fermi velocity, o~ x and o~ z are Pauli ma- 
trices, r\ is the dimensionless strength of the magnetic 
barriers, and V(x) is a non-uniform electric potential. 
A 2 x 2 unit matrix is implied in the last term. Hq is 
certainly only valid within the bulk energy gap. The 
first term is the Weyl Hamiltonian of the bare edgei 11 ! 24 
We neglect higher-order spatial derivatives, which lead 
to non-linear terms in the dispersion. An electric field 
perpendicular to the layers induces an additional Rashba 
spin-orbit-coupling term Hr = —iao~ v d x . However, this 
term can be absorbed into the Weyl term by a rotation 
in spin spaced The upper (lower) sign of the third term 
refers to parallel (antiparallel) exchange fields in the two 
barriers. The parameter rj can be written as 



V 



9^bB cxc L 

flVF 



(2) 



where g is the g- factor, \ib is the Bohr magneton, B exc 
is the exchange field, and L is the width of the magnetic 
barrier. Equation |T]) is obtained in the limit B cxc — > oo, 



L 0, keeping B CXC L constant. For the potential V(x) 
we take 



5 vy2 

V(x) = { -eV g 



for x < 0, 

for < x < W, 



(3) 



-eV/2 for x > W, 



where V is the bias voltage and V g is the gate voltage. 

Since the time-independent Schrodinger equation re- 
sulting from Hq is of first order in spatial derivatives, 
we can solve it by means of a (non-unitary) "spatial- 
evolution operator." Multiplying the Schrodinger equa- 
tion by a z from the left, we obtain 



d x ip = i ■ 



V(x) 



hvp 



a z ip(x) 



- rj <j v d(x) tjj(x) T 8{x - W) rp(x), (4) 
which for x > xq is solved by 



ip(x) 



Si- exp 



dx' 



E - V(x') 



flVF 

r/o~ y S(x')tvo- v S(x' -W) 



ip(x ), (5) 



where is a spatial-ordering directive; operators act- 
ing on ip{ x o) are ordered with their spatial coordinates 
increasing from right to left. We can obtain the boundary 
condition at the barrier at x = by making xq infinites- 
imally negative and x infinitesimally positive, 



<Mo+) 



(cosh i] — sinh ?^) ^(0 ). (6) 



Since the Schrodinger equation is of first order in spatial 
derivatives, there is only a single boundary condition. 25 
Analogously, we find the boundary condition at the other 
barrier, ip(W + ) = (cosh r\ =F a v sinh rj) ip(W~). 

To study transport, we assume the states in the leads 
to be filled up to the chemical potential fi, measured 
relative to the Weyl nodes (band crossing points), which 
are shifted by the potential ±eV/2. We take fj, to be the 
same in the two leads since they are parts of the edge of 
the same quantum well. 

In the following, we will need the eigenstates of the de- 
coupled dot, which corresponds to the limit r/ — > oo. The 
solution is straightforward: The Schrodinger equation for 
the eigenspinors ipv{x) to eigenenergies E v reads 

- ihvF<y z d x ip v - eV g tp v {x) — E y tp u (x) (7) 

with the boundary conditions 

(1 + ^)^(0)^0, (lT<r y )MW) = 0, (8) 

which give 

M0) = -i^t(O), MW) = ±%mw). (9) 

The normalized solutions are 

ex P {* ~^7^ 
-i exp 



tpu(x) 



(10) 
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with 



E v = -eV q 



■nhvp 

1/2 for parallel exchange fields, 

for antiparallel exchange fields. 



(11) 



where v can assume any integer value. Note that the 
spectrum is an equidistant ladder, which is, unlike for the 
harmonic oscillator, unbounded from both above and be- 
low. The level spacing is Eq :— irhvp/W. Of course, the 
discrete spectrum only exists inside the bulk gap and is 
only equidistant as long as non-linear terms in the disper- 
sion of the edge states can be neglected. In the absense 
of a gate voltage, the spectrum is symmetric with respect 
to zero energy for both orientations of the exchange fields 
due to the particle- hole symmetry of Hq. For antiparallel 
exchange fields, one eigenstate has energy zero. 

Finally, the Coulomb interaction is described by the 
particle-hole-symmetric term iJj nt = (e 2 /2C) (An) 2 , 
where C is the capacitance of the quantum dot and An 
is the excess number of electrons on the dot compared to 
neutrality. In terms of the number operators n u = c\,c v 
of single-particle dot states \v), we have 



E 



1 

n v -- 



(12) 



III. TRANSPORT THROUGH A 
NON-INTERACTING DOT 

Neglecting the electron-electron interaction, the cur- 
rent through the QSH quantum dot can be expressed in 
terms of its transmission coefficient T(E) by the Lan- 
dauer formula,— 



/ = 



fi+eV/2 



fi-eV/2 



dET{E). 



(13) 



Since T{E) does no depend on the bias voltage V in our 
case, the differential conductance is simply given by 



dl 
dV 



T(p + eV/2) + T(fi - eV/2) 



(14) 



The transmission coefficient T{E) is obtained from the 
transfer matrix T = 1~rT&o{Tl of the device, where the 
three factors are the transfer matrices of the right bar- 
rier, the dot region, and the left barrier, respectively. 
Since the right-moving (left-moving) electrons have spin 
up (down), Eq. ([6]) implies that the transfer matrix for 
the left barrier is 



7l = cosh n — a v sinh r\ = 



coshn i sinhn 
-i sinhn coshn 



(15) 



(The transmission coefficient of a single barrier thus 
equals 1/ cosh 2 n.) For the right barrier we get 



Tr = 



coshn ±i sinhn 
=Fi sinhn coshn 



(16) 



The transfer matrix for the dot region can be infered from 
the spatial-evolution operator in Eq. ([5]) for xq = + , 
x = W~, 



r dot = exp(i^±^^W 








(17) 



The transmission coefficient T = \t\ 2 is obtained in the 
standard way by solving 



TftTdotTL 



for t. The result is 



T(E) = 



1 



1 + trg^ 



.2 ( E+eV g 



hvf. 



W) sinh z 2n 



(18) 



(19) 



where trg = cos (sin) for parallel (antiparallel) exchange 
fields. The transmission is maximal whenever the en- 
ergy E coincides with an eigenenergy E v of the decou- 
pled dot. In this situation, the transmission coefficient is 
unity so that the device is perfectly transparent on res- 
onance. The width of the maxima is controlled by the 
strength of the magnetic barriers, n, with the width get- 
ting exponentially small for large n. This is more easily 
seen by writing the transmission coefficient as a series, 



T(E) = 



En 



2l 



2vrcosh2n ^ (E-E„) 2 +j 2 ' 



where the width of the Lorentzian peaks is 

Eq 

7 := — lncothn. 



(20) 



(21) 



An analytical expression for the current can be obtained 
by inserting T(E) into Eq. (JT3J) . We defer this to the 
next section. 

This is not required for deriving the differential con- 
ductance dl /dV, which can immediately be read off from 
Eqs. (fbfj) and (fT5|) . Note that the chemical potential /i 
and the gate voltage V g only appear in the combination 
eV g + /i since only energy differences between the dot 
and the leads enter. We can thus set [i = without 
loss of generality. dl/dV is plotted in Fig. [5] as a func- 
tion of gate and bias voltages. The figure clearly exhibits 
the suppression of the low-bias conductance off resonance 
due to Pauli blockade. The Pauli-blockade diamonds are 
filled in for strong coupling to the leads (weak magnetic 
barriers), and for n — > we recover the constant con- 
ductance dl/dV = e 2 /h of an open channel. dl/dV is 
periodic in the gate voltage with period Eq = -kHvf/W 
and in the bias voltage with period 2Eq . Also, going from 
parallel to antiparallel exchange fields has the same effect 
as shifting eV g by half a period, Eq/2. 
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FIG. 2: (Color online) Density plots of the differential con- 
ductance dl/dV through a QSH quantum dot with parallel 
exchange fields as a function of gate voltage V a and bias volt- 
age V in the absence of electron-electron interaction. The 
dimensionless strength of the magnetic barriers was set to (a) 
r\ = 1 and (b) r\ = 0.1, respectively. Light colors denote large 
dl /dV. (c) Linear conductance G = dl / dV\v=o as a function 
of the gate voltage V g . The level spacing of the decoupled 
dot, Eo = nhvF/W, is used as the energy unit. 



Additional insight that will also prove useful for the 
interacting case can be gained by writing the current in 
the form of the MW formulaM 



/ = 



fj,+eV/2 



ti-eV/2 



dE ImTrf(-2)G r et (E), 



where 



r = r lj (r L + r fl )- 1 f fl 



(22) 



(23) 



f L and f R are matrices of coupling functions to the left 
and right leads, and (^^(E) is a matrix of retarded non- 
interacting Green functions. This form is general under 
the condition that T L and T R differ at most by a scalar 
factor.— We will use the basis of decoupled dot states 
| v) . Equation (|13[) can be written in this form by taking 



GT/{E) = 



E-E v 



Vrct : 
^0 



(24) 
(25) 



where the self-energy due to the coupling to the leads is 
independent of v. Hence, we can write 



'-s 



d£r£(-2)Im— — — (26) 



with T = r L r fl /(r L + T R ), which for our case of sym- 
metric coupling implies V = T L /2 = T R /2. Comparison 
with Eqs. ([13]) and ([20) yields 



r = 



1 



^0 



Hvf 

2W cosh 2??' 

Hvf 



-i 7 



In coth?y. 



(27) 
(28) 



We note in passing that the same decay constant 7 is 
found by solving the time-dependent Schrodinger equa- 
tion for the coupled quantum dot under boundary condi- 
tions implying that no probability flows into the dot. 

In the standard MW formula the non-interacting self- 
energy is given by — i(T L + T R )/2 = — 2ir^ Our results 
contains a crucial difference: We have 



yirct 

^0 



-1 In 

2W 



2W 



2W 



(29) 



which does not coincide with — 2iT. This is because the 
MW formula was derived for a tunneling Hamiltonian, 14 
whereas our Hamiltonian Hq is not of tunneling form. 
Our model should approach a tunneling Hamiltonian in 
the limit of weak coupling, i.e., of strong magnetic barri- 
ers, 77 3> 1. In this limit we obtain 



r 



w 6 ' 



w 



(30) 



so that we indeed recover the standard relation between 
the self-energy and the coupling to the leads. 
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It is remarkable that the current can be given in a gen- 
eralized MW form, Eq. (j2"6"]l , even when the description 
in terms of a tunneling Hamiltonian is not valid. The dif- 
ference between a real-space first-quantized Hamiltonian 
and the corresponding approximate tunneling Hamilto- 
nian has been discussed by Appelbaum^ 6 . and more re- 
cently by Patton^Z, Both authors point out that the cur- 
rent calculated from the approximate tunneling Hamil- 
tonian is only correct to second order in the tunneling 
amplitudes, i.e., to first order in F j 26 i 27 We have here 
obtained a MW-type current formula for a real-space 
Hamiltonian that is exact for any coupling, at least for 
non-interacting electrons. We can also turn the argu- 
ment around: The QSH quantum dot can be described 
by a tunneling Hamiltonian with dot-lead couplings T a 
if we use the renormalized self-energy, Eq. 



IV. 



tonian 



TRANSPORT THROUGH AN 
INTERACTING DOT 



In the present section we consider the effect of electron- 
electron interaction in the form of a charging energy H ln t , 
see Eq. (fl2"|) . Since this term preserves particle-hole sym- 
metry, the current is still odd and the differential con- 
ductance is even in the bias voltage. 

To obtain the current for the interacting case, we re- 
place the non-interacting Green function G et in the gen- 
eralized MW formula by the interacting one, G rct . In 
principle, we should start from the general expression 
(f2"2"| . which allows for the Green- function matrix being 
non-diagonal in the dot single-particle basis We 
will argue below that it remains diagonal in the presence 
of Hi n t ■ A simple argument for this is that H- mt does not 
induce any scattering from one state | v) into another and 
thus leaves the self-energy diagonal in v. In this case, we 
can replace Eq. (|26|) by 



r H+eV/2 



(31) 



The simplest nontrivial approximation for G ret involves 
a Hartree-Fock decoupling of the interaction term. In 
this approximation, the spectral function A U (E) — 
-2ImG™ t (£) still is a Lorentzian at a shifted energy, 
which includes the average Coulomb interaction with the 
other electrons on the dot. 28 This can only be valid if the 
spectral function is dominated by a single peak. We now 
show that this is generally not the case. 

To go beyond the Hartree-Fock approximation, we em- 
ploy an equation-of-motion approach. The crucial as- 
sumption is that the effect of the coupling to the leads 
is not changed by the interaction. This is motivated by 
the fact that all dot single-particle states couple equally 
to the leads and that Slater determinants of these states 
remain eigenstates of the interacting dot. We implement 
this approximation by using a non-hermitian dot Hamil- 



Kdot = /XEy - i 7) ti.„ + H hlt . 



(32) 



The equation of motion for the full Green function is then 



(E - E v + i 7 ) GT V , (E) = «W + j: j dt e iEt / H 
x ({e^*/ ft [H inu c v ] e-^*/ ft ,4,}) , (33) 

where cj, {c u ) creates (annihilates) an electron in dot 
state \v). The solution of Eq. (|3"3")l is relegated to Ap- 
pendix [A] The result for the diagonal components of the 
Green function is 



...,m, y _ i ,m u -\-\ ,. . .— 

P v {. . . ,m v -i,m v+ i, . . .) 
X E-E v -% K-i)+i 7 
where the sum is over all with fi ^ v and 
P„(... ) m v _i,m y +i,-) 



- (34) 



n 



1 — Up, for m M = 



for 777, 



1 



(35) 



is the probability of the dot occupation-number state 
| . . . , m v -i, m„+i, . . .) without reference to the occupa- 
tion of single-particle state \v). Introducing the prob- 
abilities P(. . . , 777oi m ij • ■ •) °f the full dot occupation- 
number states | . . . , mo, mi, . . .), we can write 

P v {. . .,m v -\,m u +i, ■ ■ ■) 

= P(. . .,m v -x,0,m u+ i, . . .) 

+ P(...,m v -i,l,m v+ i,...). (36) 

It is also argued in Appendix[A]that the off-diagonal com- 
ponents G*j£, (E) vanish in the stationary state since they 
would require off-diagonal components of the stationary 
reduced density operator. 

The retarded Green function G^(E), Eq. ([M]). has 
a multiple-pole structure with one pole for every possi- 
ble excess electron number Aii^„ := J2^i/( m n ~ 1/2) in 
the other single- particle states ^ The energy 

of each pole contains the interaction energy of an elec- 
tron in state \v) with the excess electron number An nl/ . 
The weight of each pole is the probability to find that 
value of At7^„. This structure is highly plausible as the 
generalization of the two-pole structure for the Ander- 
son models In the non-interacting limit, e 2 /G — > 0, the 
peaks merge into a single one at the single-particle en- 
ergy. The interacting spectral function 

l 

A V {E) = ^2 P v {- ■ ■ ,m v -x,m v+ i, . . .) 

...,m 1 y_i,m 1 >-j-i,...=0 

2 7 



e-e v -%y.& v K-l) 



{37) 



+ 7 2 



6 



is a sum of Lorentzians, all of the same width 7, thanks 
to our main approximation. 



Inserting the spectral function into the generalized 
MW formula (|3"TT) we obtain the current 



J 



I = 



- / dETy AJE) = tanh — - > > 



'fi-eV/2 

x ^arctan 



P u (. . . ,m v -i,m v+ i, . . .) 



v ...,m |y _i,m iy+1 ,...— 



E„ 



E 



arctan 



C A^/v V,""M 2/ 



•it - ^ 



7 



x arctan 



IU7 



P„(. . . ,mu-i,m v +i, ■ • ■) 

v ...,m lJ -i,rnv + i,...=Q 

{^ 1/2 }-^ + #E^K-i)-M + f 



7rfe£F / 1/+1/2 



arctan 
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(38) 



r 



where the upper (lower) expression in curly braces cor- 
respond to parallel (antiparallel) exchange fields, see Eq. 
(fTTj) . Unlike in the non-interacting case, the differential 
conductance dl/dV is not simply given by the difference 
of the integrand at /1 + eV/2 and [i — eV/2, since the 
probabilities P v generally depend on the bias voltage. 

The next task is to find the probabilities P v . The sim- 
plest approximation would be to neglect correlations be- 
tween occupations of single-particle states by writing 



n 



for 

for to,, 



(39) 



Unlike the Hartree-Fock approximation, this retains the 
correct multiple-pole structure, thus transport features 
resulting from them appear in the right place in the 
conductance plots. However, neglecting correlations due 
to Coulomb repulsion overestimates charge fluctuations. 
This means that side-peaks in the spectral function carry 
too much weight. 

A better approximation is obtained by realizing that 



the probabilities P of many-particle states, which deter- 
mine the peak weights P v through Eq. (f36|) . are the sta- 
tionary solutions of a Pauli master equation, i.e., a mas- 
ter equation for the diagonal components of the reduced 
density operator. At least for a tunneling Hamiltonian, 
an exact Pauli master equation can in principle be ob- 
tained by eliminating the off-diagonal components^ - — 
In practice, one has to rely on approximations. For weak 
coupling to the leads but arbitrary interaction strengths 
one can retain only the terms of leading order in the 
coupling, which constitutes the sequential-tunneling ap- 
proximation. 

In the strict sequential-tunneling approximation, the 
broadening of dot states due to dot-lead coupling is ne- 
glected so that the probabilities show unphysical jumps 
as a function of bias voltage. We therefore include the 
broadening 7 beyond sequential tunneling^ Denoting 
the occupation-number states by to :— \ . . . , too, mi, . . .), 
the sequential-tunneling rates R^^,.^ vanish unless to 
and to' differ in exactly one of the occupation numbers 
nifi. If \v) is the corresponding single-particle state, the 
rate is 



J 



Rfh- 



R" 



2l 



rf i-eV/2 f ^+eV/2\ 

) * [E-E v -% K - I) ] 2 + 7 2 



R ( . ^{^ 1/2 }-^ + f^(r 

— 7r — arctan — 

7T V 7 



arctan 



TThvp J v+1/2 
W 



(40) 



7 



for m u = and m' u = 1 and 

\ dE 2 7 



D DOUt r> 



u-eV/2 J»+eV/2j 2* [E-E u - §J2 ^ v ( mtl -i)] + 7 



Ro / . ^ { % 1/2 } - eK, + g K 4l_-_M+j 
— 7r + arctan — 

1" \ 7 

, t ^ { V+ l' 2 }-eV 9 + b K - |) - M - f \ ^ 
+ arctan — (41) 

7 J 

I 



for m u — 1 and m' u — 0, respectively. The notation signi- 
fies that the integrals share the same integrand. The rate 
for tunneling into (out of) state \ v), R™ (R° ut ), is propor- 
tional to the spectral weight of this state below (above) 
the electrochemical potentials in the leads. Rq is a char- 
acteristic rate further discussed below. In the present 
approximation, Ro drops out of the solution for the sta- 
tionary state. The Pauli master equation now reads 

j P(rh) = J2 [Ra'^tA P{m) - Ra^A' P{m)] . (42) 

rh' 

For the stationary state the left-hand side vanishes. For 
practical calculations, the Fock space of the dot is trun- 
cated by restricting the single-particle basis to -/V sta te 
states. We have checked that increasing this cutoff does 
not lead to significant changes in the numerical results. 

Like in the non-interacting case, the chemical potential 
and the gate voltage only appear as eV g + \i in Eqs. (1551) . 
(|4"0)) . and (|4"Tj) so that the current only depends on this 
combination and we can set [i = without loss of gener- 
ality. Moreover, the current remains an odd function of 
the bias voltage and an even function of the gate voltage 
(for /i = 0). The current remains periodic in eV g but the 
period is changed to Eq + e 2 /C, since a shift of eV g by 
E + e 2 /C in Eqs. ([55]). (gDJ, and flU]) can be compen- 
sated for by shifting the ladder of single-particle states 
by one level spacing Eq , taking into account that this will 
also increase or decrease the excess charge in the station- 
ary state by one unit. Moreover, the change from parallel 
to antiparallel exchange fields is equivalent to a shift of 
eV g by half that period. For antiparallel exchange fields 
and \i = eV g — 0, a conductance peak remains pinned 
at zero bias, due to the particle-hole symmetry of the in- 
teracting model. On the other hand, there is no reason 
for the differential conductance to remain periodic in the 
bias voltage. 

Our approach is not only valid for weak coupling to 
the leads, by construction, but also becomes exact in 
the limit of weak interactions, for any coupling^ since 
then the peaks in the spectral function merge. Moreo- 
ever, the limit of strong coupling to the leads (77 — > 0) is 
also correct: The broadening of the peaks in the spectral 
function becomes 7 3> eV so that the density of states in 
the dot becomes constant and equal to 1/Eq = W/irhvF. 



'v = -3*' ^ V v 1 1 T 

v = -2 




4 



energy £ (E Q ) 



FIG. 3: (Color online) Spectral function A V (E) for sev- 
eral states \v) with single-particle energies close to the Weyl 
node. The parameters are eV = 2Eo, e 2 /C = 0.6 -Bo, and 
eV g = 0.25 -Bo, chosen so that there are no accidental degen- 
eracies. The exchange fields in the barriers are parallel. The 
level spacing Bo = ttIivf/W is used as the energy unit. The 
energies of the strongest peaks are marked by triangles. 



Thus the current becomes I = (e 2 /K)2'k(TW/'kvf)V = 
(e 2 /h) V, the result for an open channel. Note that this 
holds for arbitrary interacting strength. 

Figure |3J shows the spectral functions A V {E) for various 
v. The spectral functions for v = —3, —2, 1, 2, which cor- 
respond to large single-particle energies \E V \, have sim- 
ilar shapes and mainly differ in shifts in energy. Their 
three-peak structure shows that three charge states of 
the dot dominate. Their shape is similar since the charge 
fluctuations mostly take place in other states, namely in 
v = —1,0. The spectral functions for v = —1,0 have 
only two dominant peaks because the charge in the other 
states assumes mainly two distinct values; the fluctuat- 
ing occupation of the state \v) itself does not matter, 
see Eq. (|37|) . Also note the phase change of the pattern 
for increasing v. The spectral functions A V {E) for large 
v have an additional energy shift of e 2 /C compared to 
those for small v since for large v the total occupation of 
the other states is larger by one than for small v. 

Figure^] shows the differential conductance for various 
interaction strengths and otherwise the same parameters 
as in Fig. a). With increasing interaction strength the 
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FIG. 4: (Color online) Density plots of the differential con- 
ductance dl/dV through a QSH quantum dot with parallel 
exchange fields as a function of gate voltage V g and bias volt- 
age V. The dimensionless strength of the magnetic barriers 
is n = 1 (7 = 0.086689 Eq). The interaction strength is (a) 
e 2 /C = 0.1 So, (b) e 2 /C = 0.6 £0 (same as in Fig.[3]), and (c) 
e 2 /C = 1.2 Eo (note the shifted vertical axis). 




eV s /E 

FIG. 5: (Color online) Density plot of the differential con- 
ductance dl/dV through a QSH quantum dot with parallel 
exchange fields as a function of gate voltage V g and bias volt- 
age V for strong electron-electron interaction e 2 /C = 4.3 -Eo. 
The dimensionless strength of the magnetic barriers is n — 1. 
Other than in the previous figures the logarithm of dl /dV is 
plotted, which enhances the contrast at small dl/dV. Also 
note the extended voltage scales. 



Pauli-blockade diamonds morph into Coulomb-blockade 
diamonds and become larger, in agreement with the 
gate- voltage period Eq + e 2 /C. Also, the sequential- 
tunneling peaks defining the Coulomb-blockade thresh- 
old gain more weight; recall that in the non-interacting 
case dl/dV was periodic in the bias voltage V. 

The differential-conductance plot in Fig. |UJb) corre- 
sponds to the spectral functions A V (E) shown in Fig. [3] 
The plots are not easy to compare, since the spectral 
functions are calculated for fixed bias voltage V, but one 
can say that the splitting of the peak in A U {E) due to 
interactions is reflected by a splitting of the sequential- 
tunneling peaks in dl/dV. Importantly, weaker conduc- 
tance peaks are visible inside the Coulomb-blockade di- 
amonds. While a single dot charge dominates in this 
regime, other dot charges have non-zero probabilities, 
which are controlled by the tails of the Lorentzian peaks. 
These subdominant dot charges lead to satellite peaks in 
the differential conductance. 

Figure [5] shows dl/dV for strong interaction, e 2 /C = 
4.3 Eq. Interestingly, the conductance peaks inside the 
Coulomb diamond show an approximate periodicity with 
the non-interacting periods Eq in eV g and 2Eq in eV. 
Thus the non-interacting pattern re-emerges for strong 
interactions. Further inspection shows that there are two 
sets of lines, each with period Eq. The spectral- function 
peaks and tunneling processes corresponding to these fea- 
tures are illustrated in Fig. [6] The sequential tunneling 
lines are characterized by the state \v) involved in the 
tunneling and by the excess occupation number An-,„ of 
the other states. For strong interactions, at least one and 
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v=0. 



1,-1/2 
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FIG. 6: (a) Sketch of strong sequential-tunneling lines for 
the case of strong interaction, e 2 /C > -Bo- The lines sloping 
upward correspond to electrons tunneling out of the quantum 
dot. They are labeled by the number v of the single-particle 
state involved and by the excess electron number An n „ = 
X] M ^„( m M — 1/2) of the other states. The bold lines denote 
the edges of the Coulomb-blockade diamonds, (b) Illustration 
of the tunneling processes responsible for the lines marked in 
panel (a) for eV = and eV g tuned so that the process is 
resonant. 



typically only one peak resulting from A V (E) for fixed 
v lies inside the Coulomb-blockade diamond. Figure [B] 
shows that there are two families of lines corresponding 
to An-,,/ = 1/2 and An-, v = —1/2, respectively. (Addi- 
tional families exist but correspond to initial states with 
larger excess charge and much smaller probability.) 

We now turn to the consequences of spin-momentum 
locking. From the continuity equation for the charge, 
dtp + d x j = 0, one easily obtains the charge current in 
terms of the many-particle wave function 

eri<72...=±l I 

X ®a lt r a ..X X U X 2, ■ ■ •) CTi <F CT1(T2 ...(X1,X2, . . .). (43) 




FIG. 7: (Color online) Current as a function of gate voltage 
for bias voltage eV = -Bo, interaction strength e 2 /C = 4.3 Bo, 
and two strengths r\ of the magnetic barriers, which have par- 
allel exchange fields. The current is proportional to the uni- 
form z-component of the spin density. 



This expression is clearly proportional to the z-compo 
nent of the spin density p z , 

2evF 



■p*(x). 



(44) 



In the stationary state, the current j(x) = I is uniform 
so that the z-component of the spin density is also uni- 
form and given by p z = —%I/2evp. Thus the quantum- 
dot structure permits to control the magnetization of 
the whole edge with a gate voltage applied locally to 
the dot. We do not address the dynamics in this pa- 
per, i.e., we do not study how the edge reacts to a non- 
adiabatic change of the gate voltage. Figure [7] shows 
the uniform current (spin-z density) as a function of the 
gate voltage for strong interaction, e 2 /C = 4.3 Eq, corre- 
sponding to Fig. [5] For stronger magnetic barriers, the 
current and the spin density are reduced but the con- 
trast between sequential-tunneling ( "on" ) and Coulomb- 
blockade ("off") regimes is higher. 

Equation (|44j) begs the question whether the z-compo- 
nent of the spin current is related to the charge density. 
Indeed, from the spin continuity equation d t p z + d x j z = 
one obtains the spin current 

Hvf 



} dxidx 2 ... y^(x - Xi) 

ai(72...— ±1 i 
X <r 2 . ..(£1)^2, • • ■)^a 1 a 2 ...{x U X 2 , ■ ■ ■) 



2e 



p(x). 



(45) 



This result does not hold in the magnetic barriers, since 
spin is not conserved there. For the same reason, the 
spin-z current and thus the charge density are generally 
not continuous through the barriers. An exact result can 
be obtained at eV g — for both parallel and antiparallel 
exchange fields: Here we have p(x) = and thus j z (x) — 
within the dot for arbitrary bias voltage due to particle- 
hole symmetry. 
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The approximations involved in describing the cou- 
pling to the leads by a non-hermitian Hamiltonian and 
in truncating the equations of motion for the Green func- 
tions as shown in Appendix [A] may appear to be unsys- 
tematic. Insight into their meaning can be gained from 
an alternative derivation of the same final expression for 
the current. Since we calculate the stationary proba- 
bilities P(m) in an effective tunneling picture with bare 
tunneling rate Rq for each lead, it is natural to write 
down the current in the same picture. In the geneal- 
ized sequential-tunneling approximation including level 
broadening, the current through lead a is^Sr— 



-a- I > 



(46) 

where a = L = l(a — R = —1) for the left (right) 
lead. RZ^ m , denotes the part of the rates in Eqs. (|4"0")) 
and (|4ip due to lead a. For the symmetrized current we 
obtain, using an obvious notation, 



/ = 



2h ? ^ 

m v 



R,, 
-R 



for m v = 



L.out 



^fl.out for 



1 



P(m). (47) 



We reobtain Eq. for the current if we assume that 
the scale of the sequential-tunneling rates is the coupling 
function Rq = T L = T . This is plausible but was not 
a-priori obvious since our Hamiltonian is not of tunnel- 
ing form. The approximations used in the first approach 
to find the mapping of the probabilites P(m) to the cur- 
rent are thus equivalent to a sequential-tunneling approx- 
imation with additional level broadening. The first ap- 
proach has the advantages of also providing the spectral 
functions A v (E) and of being suitable for approximations 
other than perturbation theory for weak coupling. 



V. SUMMARY 



leads to an equidistant ladder of dot resonances and to 
a double periodicity of the differential conductance in 
both gate and bias voltage. The results can be ana- 
lyzed within the framework of Meir-Wingreen theory for 
arbitrary coupling between dot and leads although the 
Hamiltonian is not of tunneling form. 

In the interacting case, an equation-of-motion ap- 
proach for the local retarded Green function yields a 
multiple-peak structure in the spectral function, where 
the peaks result from distinct charge states of the dot. 
Their probabilities are calculated within a rate-equation 
approach, which is valid for high magnetic tunneling bar- 
riers. The approach also becomes exact in the two lim- 
its of vanishing tunneling barriers and of weak interac- 
tions. We have calculated and analyzed the differen- 
tial conductance from weak to strong Coulomb interac- 
tion on the dot. For increasing interaction strength, the 
Coulomb-blockade diamonds grow in size. The periodic- 
ity in the gate voltage is retained, while the periodicity 
in the bias voltage is destroyed. For stronger interac- 
tions, weak tunneling peaks appear inside the Coulomb- 
blockade diamonds. For very stong interactions com- 
pared to the spacing of the dot resonances, two copies 
of the non- interacting periodic pattern re-emerge inside 
the diamonds. 

Spin-momentum locking leads to an obvious pro- 
portionality between the charge current and the z- 
component of the spin density or magnetization in the 
whole QSH edge. This should make it possible to control 
the magnetization by a locally applied gate voltage. 
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We have considered a quantum dot in a QSH edge real- 
ized by two thin magnetic tunneling barriers with parallel 
or antiparallel exchange fields. For vanishing electron- 
electron interactions, the linear dispersion of edge states 



Appendix A: Evaluation of the equation of motion 



In the following solution of the equation of motion ([33 
we set h = 1. Equation (|33l) contains the commutator 
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2 U ° v 



-E 
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Noting that n M and R\ Yi t commute with Kdot, it follows that 



iK dot t 
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ex p(-^E (v-^W 



Next, we have to evaluate the anticommutator {e tKdott [Hi nt ,c„]e lKdatt ,c\,}. For the diagonal case v 
immediately get 
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while for v 7^ z/ we find 
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The equation of motion for the diagonal components of the Green function thus reads 



f°° B I 
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2C 



= - dti{E-E v + i 7 ) e i(«-^)*-7t ^ JJ 
Hence, the diagonal Green function is 
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The sum is over all m M with /i 7^ v. 

For the off-diagonal components, v 7^ f', inserting Eq. (|A4|) into the equation of motion (|33|) gives 



/•OO 

(E - E v + i 7) Gj£ (£) = 2* / dt e" 

Jo 



dt I 2C 
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eXP ' * 2C * ' ^ ~ + ° XP ' ~* ^ 1 
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dt (E - E v + i 7 ) e *(E-^)*-7* sin ( |_ t 
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Performing the integral, we obtain the off-diagonal Green function 

cm = t (n{ 1 ; ,, £;:!}'U 



...,mo,mi,...= \ii^v,v' 



' ' (A8) 



where the sum is now over all with p, ^ v,v' . The off- 
diagonal components of the Green function thus contain 
averages involving the off-diagonal operator c* u ,c v . In the 
stationary state such averages vanish: For the average 
in Eq. (|A8|) to be non-zero, the stationary-state reduced 
density operator p would have to be non-diagonal in the 
occupation-number basis. Specifically, superpositions of 
states different by a single electron being transferred from 



\v) to \v') would have to occur. From the decoupled dot 
Hamiltonian, such off-diagonal components of p obtain a 
rotating phase factor with frequency uj = \E u i — E u \/h =/= 
0. For the stationary state, this time dependence would 
have to be canceled by the time dependence due to the 
coupling to the leads. That this cannot happen follows 
for example from the fact that the leads do not contain 
the energy scale E = itHvf/W. 
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